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Abstract 


We define the generalized basic hypergeometric polynomial of degree N as follows: 
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Here N is an arbitrary positive integer, r and s are arbitrary nonnegative integers, the r + s parame¬ 
ters Uj and Pf. are arbitrary (“generic”, possibly complex) numbers, {a\q)^ is the g-Pochhammer sym¬ 
bol and r-i-i</>s (q^Oj cn,---) ttr;/3i,...,/Sgl g; z) is the generalized basic hypergeometric function. In this paper 
we obtain a set of N nonlinear algebraic equations satisfied by the N zeros = Cn (Sj/3; 9!-^) of this 
polynomial. We moreover identify an (TV x A^)-matrix M = M (a, P; (; q; JV) featuring the JV eigenvalues 

r 

Pn = j-^-n — 1) p[ (ofj g^“" ~ 1) : TL = 1,2, ..., N . These N eigenvalues depend only on the 

i=i 

r parameters aj (besides g and N), implying that the (TV x A)-matrix M is isospectral for variations of the 
s parameters Pj.; and they clearly are rational numbers if g and the r parameters aj are themselves rational 
numbers: a nontrivial Diophantine property. 


Keywords: basic hypergeometric polynomials, Diophantine properties, isospectral matrices, special func¬ 
tions. 


1 Introduction 


We define the generalized basic hypergeometric polynomial of degree N as follows: 
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Here N is an arbitrary positive integer, r and s are arbitrary nonnegative integers, the r + s parameters aj 
and pf. are arbitrary (generic, possibly complex) numbers, (7; g),^ is the g-Pochhammer symbol, 

(7;9)o = l> (7;9)m = (l-7) (I- 79) (1-79b • • • (1-79™“b for m = I, 2,3,... , (2) 

and r+i4>s (q^Oj CKr;/?!, ■•■,/3s; 9;-z) is the generalized basic hypergeometric function (see ([HI) below for its 
definition, of course consistent with ©; and, for instance, [I] for its properties). We actually prefer to work 
hereafter with the monic version of this polynomial, pn (a; /3; g; z) , defined of course as follows: 
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which is clearly characterized by the property 
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In this paper we report properties satisfied by the N zeros = Cn of this polynomial pn {z) 

(or, equivalently, (2:)), namely by the N roots Cn defined (up to permutations) by the formula 

PAT (a; g; C„) = 0 , n = 1, 2,..., TV , (4a) 

SO that 

N 

PN {gr,^,q]z) = '^[z - C,^[g:,^,q\N)] . (4b) 
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Notation 1.1. Let us confirm that throughout this paper the parameters r, s and N are 3 arbitrary 
nonnegative integers (except when their values are explicitly assigned, see for instance the subsections at the 
end of next section). Above (and hereafter) the short-hand notations a, /3 respectively C denote the unordered 
sets of the r, s respectively N numbers aj, /3j. respectively The N zeros are of course functions of the 
r -I- s -I- 2 parameters aj, q ^ 1 and N, i. e. C = C (n; /?; 9! -^)■ Note that we occasionally omit to indicate 
explicitly the dependence on some, or on all, of these parameters (including, systematically, the dependence on 
r and s, the values of which are considered fixed throughout this paper). Generally (unless otherwise specified) 
the indices n, to, £ run over the positive integers from 1 to TV, while the index j runs from 1 to r and the 
index k from 1 to s. Upper-case underlined letters denote TV x TV matrices: for instance the matrix M has 
the TV^ elements Mnm- We moreover assume the parameters aj, q ^ 1 to have generic (possibly complex) 
values, and the TV zeros Cn to be all different among themselves, Cn 7^ Cm for n ^ m, which is of course the 
generic case; on the understanding that otherwise our formulas might have to be understood cum grano salis, 
via appropriate limiting processes. □ 

In this paper we firstly obtain a set of algebraic equations satisfied by the TV zeros Cn of the generalized 
basic hypergeometric polynomial of order TV, additional to the algebraic equation defining these zeros, see (I4ap . 
This result—reported below as Proposition 2.1 because we were unable to hnd any mention of it in the 
literature—does not seem to us to be particularly significant, but it is instrumental to prove what we consider 
our main finding: the identification of an (TV x TV)-matrix M that features the TV eigenvalues 


Pn = (g-” - 1) n - 1) ’ ^ 2- -- ^ ° (5) 

This matrix M is explicitly defined below in terms of the r-|-s-|-2 parameters Oj, /T^., g 7^ I and TV characterizing 
the generalized basic hypergeometric polynomial pn [z), and of its TV zeros Cn (which themselves depend of 
course on the r -|- s -|- 2 parameters Uj, Pf., q ^ 1 and TV). While the TV eigenvalues clearly depend only on 
the r -|- 2 parameters 07, g 7^ 1 and TV, implying that the (TV x TV)-matrix M_ is isospectral for variations of 
the s parameters P/.', and all these eigenvalues are clearly rational numbers if g and the r parameters aj are 
themselves rational numbers: a nontrivial Diophantine property. 

The hndings outlined above are detailed in the following Section 2, and our main hnding is proven in 
Section 3. A terse Section 4 (“Outlook”) outlines possible future developments. The definitions and some 
standard properties of the generalized basic hypergeometric polynomials are reported in the Appendix, for 
the convenience of the reader and also to specify our notation; the reader is advised to glance through this 
Appendix before reading the next section, and then to return to it whenever appropriate. 

For a terse review of analogous results for generalized hypergeometric polynomials the interested reader is 
referred to [2]; of course the results reported in this reference [2] can be obtained (up to obvious notational 
changes) from those reported in the present paper by taking appropriately the g —>■ 1 limit. For results analogous 
to those obtained here, but for polynomials belonging to the Askey and g-Askey schemes, see i i; and let 
us recall in this connection that the polynomials belonging to the Askey, respectively the g-Askey, schemes 
(see for instance 0) are also special cases of generalized hypergeometric, respectively basic hypergeometric, 
functions, with indices r < 3 and s < 3, in which the polynomial variable is however related to parameters of 
the hypergeometric function rather than to its argument] hence the results reported in 0 a are not special 
cases of those obtained and reported in [2] and below. 

Finally, let us mention that the properties of the zeros of polynomials are a core problem of mathematics to 
which, over time, an immense number of investigations have been devoted. Nevertheless new hndings in this 
area continue to emerge, see, for instance 0 [7] [5] [5] [TU] [H] (and of course 0 a 0 )- 
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2 Results 


The generalized basic hypergeometric function r+i4>s (ooj ...,l3g;q;z) is defined as follows (see for 

instance eq. (1.2.22) in [I], or eq. (0.4.2) in [5]): 
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where the g-Pochhammer symbol (a; q)^ is defined above, see ([U, and the rest of the notation is, we trust, self- 
evident. Clearly if one of the r-|-1 parameters aj coincides with a negative integer power of g, say (without loss 
of generality) ap = q~^, and all the other r -I- s parameters aj and have generic (possibly complex) values, 
the series in the right-hand side of the definition (|6]) of the generalized hypergeometric function terminates at 
j = N (since clearly (g“^; g)^ = 0 for p = TV -|- 1, -I- 2,..., see ([2])). This leads to the definition ([T]) of the 

generalized basic hypergeometric polynomial Pn (z) (of degree N in z; below, we rather work with the monic 
version pn {z) of this polynomial, see ©)■ 

The first result of this paper—proven in the Appendix—is the following 

Proposition 2.1. The N zeros = C,j^{ar,P]q]N) of the generalized basic hypergeometric (monic) 
polynomial pn (o; P', g; z) , see (P), dSil) and (f4a|l . satisfy the following set of N algebraic equations: 
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Here the s quantities bj. (/3) , respectively the r quantities aj (a), are defined, in terms of the s parameters pj^ 
respectively the r parameters Uj, by 621) respectively 63 • Note that the formula © becomes a bit simpler 
in the “balanced” case (i. e., if r = s). Also note that, if instead s > r, some terms in the sum over the index 
j have a vanishing value (for j = s — roTj = s — r — 1). □ 

Our main finding—proven in Section 3—is the following 

Proposition 2.2. Let the {N x A^)-matrix M be defined, componentwise, as follows in terms of the r-|-s-|-2 
parameters aj, P^, g, N characterizing the generalized basic hypergeometric polynomial and of its N zeros 
C„ = Cn (ffi TV), see ©, dSil) and 6il): 

Mnn = (-1)^1 (<7 - 1)"<?«(1,C) + - l)V(fc + 1,C) - {q^ - l)V(ft,C)] I 
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+ E + s + 1 - r, 0 - - l)^9n{j + s - r, C)] 

T=i 


+ (-l)q g-^(g«-”+i - l)/„(s - r + 1, C) - {q^-^ - l)Us - r, Q 


-f g ^(g-?+®+i - l)/„(2 + s + 1 - r,C) - (g^^'* - l)/n(j + S - r,C) 

T=i 

n = 1,2,...,TV, 


(8a) 


3 


















{q - l)^/nm(l,C) 


Mnn^ = (- 1 ) 


s+1 




(C -C 

\Sn Sm/ 


+ E - Iffnmik + 1,0 - - l)Vn-(fc, 0] 


fc=l 


+(-ir 


C: 


(C -C 

vSn Sm/ 


<7 (0 - 1) /nm(s - r + 1,C) - ( 9 ® - l)^/„m(s -r,C) 


+ - l)Vn,n(j + S + 1 - r,0 - (0+^-^ - ifM + S - r,0] , 
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n,m = 1,2,..., N, n^m, 


(8b) 


where the quantities fnmip, C) gn{p, C) ^'^e defined by (1401) . Then the N eigenvalues /j,„ of this matrix are 
given by the following neat formulas: 


i=i 


(ajq^-^-l), n = l,2,...,N .n 
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The following corollary, which is an immediate consequence of Proposition 2.2, yields, via the definition 
([5]) of the {N X A^)-matrix M and the expression ([9]) of its N eigenvalues /i„, a number of algebraic formulas 
satisfied by the N zeros of the generalized basic hypergeometric polynomial of order N, see o or ( 1 ^ and 
((iai) . 

Corollary 2.1. 

N 

Trace [(M)n = E , P = 1,2,3,..., (10a) 
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N 

Det (M) = n ° ( 10 b) 
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Remark 2.1. For fixed q and N, the N eigenvalues /i„ of the {N x A^)-matrix M (see (| 8 ])) depend only 
on the r parameters aj (see (O), while the matrix M_ depends on the s + r parameters /3j, and aj via the 
dependence of the parameters bk respectively aj on respectively on aj (see (SZl) respectively ®) and via 
the dependence of the N zeros on the parameters /3^ and aj. Hence the {N x lV)-matrix M is isospectral 
for variations of the s parameters P).. And note moreover that the N eigenvalues are rational numbers if 
the r parameters aj are themselves rational numbers, and q is also a rational number: a nontrivial Diophantine 
property of the {N x iV)-matrix M. □ 

Remark 2.2. All the above results are of course true as written only provided the N zeros Cn 
different among themselves; but they remain valid by taking appropriate limits whenever this restriction does 
not hold. □ 

Remark 2.3. Immediate generalizations—whose explicit formulations can be left to the interested reader— 
of Propositions 2.1 and 2.2 obtain from these two propositions via the special assignment Of+p = Pg+p 
p = 1, ...,u with u an arbitrary nonnegative integer such that both f = r — u and s = s — u are positive 
integers. These propositions refer then to the N zeros of the generalized basic hypergeometric polynomial 
pNiai,...,af;Pi,...,l3g;q;z) = Pn (ai,... ,af,.. ■ ,ar; , Pg,..., l3^;q; z) which depend—additionally 

to N and q —only on the f+s = r+s —2u parameters aj with j = 1,..., f = r—u and /3j. with fc = 1,..., s = s—u, 
but feature quantities bk and aj (see dni) and ([49|)) that depend on the 2 + r + s parameters N, q, aj with 
j = 1 ,..., r and Pk with fc = 1 ,..., s. □ 

Let us end this Section 2 by displaying explicitly the above results for small values of the integers r, s and 
u (see Remark 2.3). 


2.1 The case r = s = l, u = 0 

If r = s = 1 and u = 0, the generalized basic hypergeometric polynomial o is given by 
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Let C„, where n = 1, 2,..., be the zeros of P/v(ai; /3i; <z; z). Then, by Proposition 2.1 
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+ (-^+aiCn'7 (Cn9^ “ Cm) = 0 ) n = 1, 2,..., iV . 


Define an (TV x A^)-matrix M componentwise as follows: 

Mnn = {q-lfgn{l,0 - 1 - y+C„('?"^+ ai) + C) ^ 
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, n,m = 1,2,..., N, n ^ m 


(13b) 


where the quantities fnm{pX) and g„(p, C) are defined by (jiDl) . By Proposition 2.2, the N eigenvalues /r„ of 
this matrix are given by 

Mn = -(9“"-l)(ai9'^“"-l), n = l,2,...,iV. (14) 

By Corollary 2.1, the trace of the matrix M is given by 


Tr(M) = - q-^^-^) + i±Xll - (1 _ _ jy - 1 

0 ^ — 1 a — 1 


(15a) 


and the determinant of M 
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det(M) = (-1)^ n [(9“” - l)(«i9'^“” - 1) 
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(15b) 


2.2 The case r = 2, s = l, u = 0 

If r = 2, s = 1 and u = 0, the generalized basic hypergeometric polynomial o is given by 
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Let C„, where n = 1, 2,..., A^, be the zeros of P/v(ai, 02 ; Pi] 9 ; z) and let oi = oi (a) = oi + 02 , 02 = 02 (o) = 
0102 (see (14^ 1. Then, by Proposition 2.1, 
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Moreover, by Proposition 2.2, the N x N matrix M defined componentwise by 
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has the eigenvalues 

Mn = -9 

By Corollary 2.1, the trace of the matrix M is given by 

^-N 
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and the determinant of M 


N 

det(M) = (-1)"^ n [9”'^+"(9“" - l)(ai9'^-" - l)(a29'^-” - 1) 
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2.3 The case r = s = 2, u = 0 

If r = s = 2 and u = 0, the generalized basic hypergeometric polynomial ([T]) is given by 
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Let Crt; where n = 1,2,...,A^, be the zeros of PN(oii,a2\ Hi, let oi = ai (a) = oi + 02 , 02 = 

02 (a) = Q!rQ! 2 , (see dMI), hi = bi {P) = /?! + 132, &2 = b 2 (^) = /3 i^2 (see (1471) 1. Then, by Proposition 2.1, 
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Moreover, by Proposition 2.2, the N x N matrix M defined component-wise by 
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has the eigenvalues 


(24) 


/i„ = -( 9 -” - - l)(a 2 g'^-” - 1), n = 1, 2,..., TV. 

By Corollary 2.1, the trace of the matrix M is given by 
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Tr(M) = - ^ [(q-" - l)(aig^-" - l)(a 2 (z'^-” - 1) 


(25a) 


and the determinant of M 

det(M) = (-1)^ n [(«■” - - 1) 


N 


(25b) 


2.4 The case r = s = 2, m = 1 

If r = s = 2 and u = 1, then P 2 = 02 and the generalized basic hypergeometric polynomial ([T]) is given by 
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Let Cn, where n = 1, 2,..., be the zeros of PAr(ai, a 2 ; Pp, q; z) and let m = ai + 02 , 02 = 01102 , (see (Id^ l. 
bi = Pi + a 2 , 62 = /3 i02 (see (liTI) '). Then, by Proposition 2.1 and Corollary 2.1, the zeros ( satisfy algebraic 
relations (|2^ and the matrix M defined in terms of the zeros ( by formulas (??) has the eigenvalues 


the trace 


and the determinant 


Pu = -(«■” - l)(oig'^-” - l)( 02 g'^-” - 1), n = 1,2,..., iV, 
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Tr(M) = - E [( 9 ”" - - 1 ) 
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det(M) = (-1)^ n [(q-" - l)(oig^-" - l)(o 2 (?^-" - 1) 


n—1 


(27) 
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3 Proof of Proposition 2.2 

Let the t-dependent monic polynomial, of degree N in z, be characterized by its N zeros Zn (t), 
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i’Niz.t) = n [z-Zn{t)] , 


n—1 


and by the N coefficients Cm (t) of its expansion in powers of z, 
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tpN {Z,t) = + E [cm (t) Z^ "*] . 


(29a) 


(29b) 
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It is plain that these two representations, (I29al) and (I29bl) . are consistent, implying a uniquely defined expression 
of the N coefficients Cm (t) in terms of the N zeros z„ (t) , and an expression—unique up to permutations and 
of course explicitly known only for N < 4—of the N zeros (t) in terms of the N coefficients Cm (t) ■ 

The next step is to assume that this t-dependent polynomial ijjj^ (z, t) satisfy the following linear Differential 
g-Difference Equation (DgDE): 


dt 


If (^“i) 
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If i^dk/q) 
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(30) 
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where the operator A^, acting on functions of the variable z, is defined as follows: 


= (7 (5 - 1) with 5 f{z)= f {q z) 


(31a) 


(see (j44bp i. Note that this implies 

A^ zP = {jqP -1) zP . (31b) 

Let us first of all check whether this is consistent with the fact that '(/'at {z, t) is a polynomial of degree N 
in z. Since the operator A-y, when applied to a power of z, does not change that power but only multiplies it 
by a number, see (|31bp . it is plain that this is guaranteed by the fact that A^-n = (1 — 1) z^ = 0 (see 
(j31b(l l: while of course clearly Aic = 0 (where c indicates any constant, i. e. a quantity independent of the 
variable z). 

Next, let us see what the DgDE ((5D1) implies for the coefficients Cm (t), see (|29bl) . Via (I29bl) and (|31bl) it 
clearly amounts to the following autonomous linear system of ODEs: 


(g^—+1-1) X{{P^q^-^-l) 

k^l 

r 

^s-r){N-m) - l) - l) 


Cm —1 (0 
Cm (0 




m = 1 , 2 ,..., A , with cq = 1 . 


(32a) 


Of course here and below a superimposed dot denotes a t-differentiation. 

Note the formal consistency of this system of evolution equations with the assignment cq = 1, hence with 
the fact that (z, t) is a monic polynomial of degree N in z, see (I29bl) . 

It is plain that the solution of this system reads 
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w = E [■ 


Vn exp (/r„ t ) 
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(32b) 


where the A-vectors respectively the A numbers are the A eigenvectors, respectively the correspond¬ 

ing A eigenvalues, of the (A x A)-matrix C_ with elements 

r 

Cnn = (g-" - 1) n (“J - 1) ’ ^ = 1- 2,..., A , (32c) 

1=1 

S 

Cn,n-1 = - 1) n 9 "^”” - 1) > = 2, N , (32d) 

k^l 

and all other elements vanishing: 

c = IJ-n ; (32e) 

while the A (t-independent) numbers rj^ can be arbitrarily assigned—getting thereby the general solution of 
system (I32all —or can be adjusted to fit the initial data Cm (0) so that 
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Vn K 


TO = 1, ..., A 
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—getting thereby the solution of the initial-value problem for system (I32all . The triangular character of the 
(A X A)-matrix C_, see (|32cl) and (j32dp . implies that its A eigenvalues coincide with its diagonal elements: 


/!„ = (g-" - 1) n 9 "^”" - 1) > ^ = 1’ 2, N . (32g) 

1=1 


Our next task is to discuss the t-evolution of the A zeros z„ (t) of tp {z, t ), see (I38al) . implied by the DgDE 































The first observation is that the equilibrium—i. e., t-independent—solution tp {^) of DgDE (l30l) is the 
generalized basic hypergeometric polynomial ppf (z) = pn{,cx', Q] z), see (I5a1) . 

i’iz) =Pn{z) ■ (33) 

This is implied by the g-difference equation satisfied by pjv {z), see (HSl) . which clearly implies that the right- 
hand side of the DgDE (1511 vanishes for ip {z, t) = Pp (z) = pjq (z) (and note the consistency implied by the fact 
that '0 (z) and piq (z) are both monic polynomials of degree N). Hence the equilibrium—i. e., t-independent— 
configuration of the N zeros z„ (t) of ip (z, t), see (I38al) . is 

Zn{t) = ZniO) = Zn = , n = 1,2, N , (34) 

where the N numbers are the N zeros of the (monic) generalized basic hypergeometric polynomial pn (z), 
see (|4a)) . 

Next, let us reformulate DgDE (1501) as follows: 

(^) '^N {z,t) = RHS {z,t) , (35a) 


RHS {z,t) = (-l)®z 40^(gz,t) - ’<pj,!{z,t) + ^6fc(^) ^ [tpj,f{q'^+^z,t) - tpj,f{q'^z,t)] 

[ k=l ^ 

r 

-b^aj(a)(-l)^ [q-^z,t) - ipN{q^+^-^z,t)] 
i=i 

by repeating, on the right-hand side of DgDE (1501) . the same development that led, in the Appendix, from (I44al) 
to (|50c|) : this of course implies that the quantities 5^ (/3) respectively Uj (a) are, here and hereafter, defined as 
in the Appendix in terms of the parameters (3^. respectively aj, see (!T7|) and dm). 

Next, let us insert in this DgDE, (I35L the representation (I29al) of the monic polynomial ip {z,t) via its zeros 
Zn (t). The left-hand side of this DgDE then reads (by logarithmic t-differentiation of (I29al) l 


(35b) 


) ^N{z,t) 


N 

-^N{z,t) 

m—1 


Zra jt) 

Z - Zm (t) 


= - ^\zm{t) 


n (t)] 


(36a) 


implying, for z = z„ (t). 


d t 


ipN iz:,t) 


N 


= -Zn (t) [Zn (t) - Zi (t)] . 

2=Zn{t) 


While the right-hand side (I35bl) of (I35al) clearly reads 


N 


N 


RHS {z,t)= (-i)®z n “ n 


. 1=1 




i-if 


k^l 


N 


N 


[g'^+iz - Zi{t)] - [q’^z - Zi{t)] 


U =1 




N 


N 


£^1 

N N 

q-N [ql+s+i-r^ _ _ Y[[q^+^-^z - z,(t))] 


3 


-H^Oj(a)(-l) 
0 = 1 


1=1 


(36b) 


(37a) 
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implying, for z = Zn (t) 


N 


RHS (z„ = (-l)®z„(t) 0 p [qZn{t) - Zi{t)] 


+y] {p) 


{-If 


k=l 


N 


N 


n - zft)] - P [f^Znit) - zft)] 


e=i 


1=1 


N 


N 


-(-1)’' q-^ n - n w~''^n{t) - zft)] 


1=1 


i=i 


Y.^Aa){-iy 

i=i 


N 


N 


g ^ n ''Zn{t) - zft)] - P - zft)] 


1=1 


(37b) 


It is thus seen that the equations of motion characterizing the t-evolution of the N zeros {t) of iji (z, t) 
read as follows (of course below a superimposed dot denotes a t-differentiation, and we omit for notational 
simplicity to display the t-dependence of the zeros): 

i„ = (-l)*+i| {q - l)fn{l,z) + ^ bk{p)t}^ [(y+i _ l)f^{k + l,z) - {q'^ - l)U{k,z) 

I fe=l ^ 

+ (-l)’'zJ - l)/„(s - r + l,z) - (g*”’' - l)/„(s - r,z) 


+ y]aj(a)(-l)^ q ^(g-J+"+i ^ - l)f^{j + s + 1 - r, z) - {q^+^ - l)/n(j + s - r, z) L (38a) 

i=i J 


where 


fn{p,z) = fn{p,Zi,. .. ,Zn) = J| — —'\,n = l,2,...,N. 


(38b) 


This is an interesting dynamical system, a complete investigation of which is beyond the scope of the present 
paper. But before proceeding with our task, let us pause and recall that the first idea to relate the zeros of 
polynomials to the equilibria of a dynamical system goes back to Stieltjes and Szego [12) . was resuscitated in 
m to identify “solvable” many-body problems (see also the extended treatment of this approach in [2]), and 
then extensively used to obtain results concerning the zeros of the classical polynomials and of Bessel functions, 
see the paper |15| where several such findings are derived and reviewed. For more recent developments along 
somewhat analogous lines see, for instance, m, m, [m, m, m- 

Here we need to focus only on the behavior of this dynamical system, (j38a|l . in the immediate neighborhood 
of its equilibrium configuration z„ (t) = = Cn) see ([34]). To this end we set 


f) = C„ + e (0 : n=l, .., N , 

implying 

in (t) = e Cn (0 : ^ = 1, N , 

with e infinitesimal. To aid the task of linearizing system (I38al) . we introduce the quantities 

N / 

qP Zn - Zl 


fnm{p,z) = fnm{p,Zl,...,ZN) = H f 

0 1 B-L \ J 


(implying fnm {p,f) = fn{p,z) [{z„ - Zm) / {q^ Zn - Zm)\, See (l38bD ), and 

N p 

gn{p,z) = gn{p,Zl,...,ZN) = fnk{p,z) 


Zk 


and we note that 


k—l.k^n ^ 


y^{p,z) = (1 - qP) gnip,z) 


{Zn-Zkf\ 


(39a) 

(39b) 

(40a) 

(40b) 

(40c) 
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and 


^{P,Z) = - 1 ) fnm{p,z) , (40d) 

\Z^ Zjji j 

where n,m = 1,2,..., N and n ^ m. 

The insertion of ansatz dMI in the equations of motion (I38ap is, to order £° = 1, clearly consistent via 
Proposition 2.1. To order e, it yields the linearized system of ODEs 


N 

i = M^-, , n = l,...,N, (41) 

m—1 

with the {N x 7V)-matrix M dehned, componentwise, by (|S]) (the diligent reader might wish to check, using 
the notation and the formulas in (|38bl) and (I40p . the relevant computation, which is standard but somewhat 
cumbersome—although certainly doable without the help of a computer). Here and hereafter the notation 
^ = £_{t) denotes the iV-vector of components (<). As for the implications to orders with p = 2,3 ,... 

of the insertion of the ansatz (IMl) in the equations of motion (I38al) . they consist of additional systems of 
algebraic equations satisfied by the zeros of the generalized basic hypergeometric polynomial Pn (z) or 
Pn {z), see ([T]), (I5al) and (Hal) , the explicit display of which is left to the interested reader (to get them it might 
be expedient to take advantage of appropriate computer packages such as Mathematica or Maple). 

The general solution of the system of linear ODEs (IM]) reads of course as follows: 


N 

m—1 


(42) 


where the N (t-independent) parameters can be arbitrarily assigned (or adjusted to satisfy the N ini¬ 
tial conditions (0)), while the numbers respectively the A'-vectors are clearly the N eigenvalues 
respectively the N eigenvectors of the matrix M . 

m = l,...,N. (43) 

But the behavior of the dynamical system dSD in the immediate vicinity of its equilibria cannot differ from its 
general behavior, which is characterized by the N exponentials exp t), as implied by the relation between 
the N zeros Zn (t) and the coefficients Cm (t) of the monic polynomial (of degree N in z) i/'at {z,t), see (H5|) . 
and by the explicit formula, see (1321) . detailing the t-evolution of the N coefficients Cm (t) . Hence the (set of) 
eigenvalues of the matrix M_, see ®, must coincide with the (set of) eigenvalues p^, see ( |32g[ ), of the 
matrix C, see (1321) . Proposition 2.2 is thereby proven. 


4 Outlook 

A follow-up to the present paper might explore the dynamical system (I38al) . as well as other dynamical systems 
obtained by analogous methods, also connected with hypergeometric or basic hypergeometric polynomials, but 
featuring more interesting equations of motion, for instance of Newtonian type—allowing their interpretation 
as “many-body models”. 

Another possible direction of further investigation might try and extend the approach and findings, reported 
in this paper for the N zeros of basic hypergeometric polynomials of order N, to the, generally infinite, set of 
zeros of (nonpolynomial) generalized basic hypergeometric functions. 
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6 Appendix A: Properties of the generalized basic hypergeometric 
polynomials 


In this appendix we report some properties of the generalized basic hypergeometric polynomial (z) = 
Pn {a;/3;q;z) (see ([T]) and (153 and Notation 1.1); and we prove Proposition 2.1. It is understood that, 
throughout this Appendix, the r + s + 2 parameters aj, /3j,, q ^ I and N have fixed values; the indication of 
the dependence upon them is omitted whenever this simplifies the notation at no cost in terms of clarity. 

Our starting point is the following g-difference equation satisfied by the generalized basic hypergeometric 
function r+i'Ps (q^Oj cti,---; (see Exercise 1.31 on page 27 of [I], and note the replacement of r 

with r + 1 and of Oq with q~^)'. 




-2 Aq-N 


n 


/q> 


.k^l 


n (*.,) 


r+i'/'s {q~^, ar] Pi,/3/,q; zq^-^) = 0 


(44a) 


_i=i 

where the difference operators A.^ and 5 operate as follows on functions / {z) of the variable z: 

A 7 / {z) = 7/ {zq) - f{z) , A^=-f5-l , 5 f{z) = f {zq) . (44b) 

Note that the operators A-^^^ and A^^ commute, 

Ai A 2 / (^) = 7 i 72 / (q^z) - (7i + 72) / {qz) + f{z) = A^^A^^ f (z). (44c) 

Via ([1]) and (15a)) this implies that the generalized basic hypergeometric polynomial of order N satisfies the 
following g-difference equation: 


Ai 


n 


/g) 


.k=l 


PN {a;^,q;z) - z A^-n 


i=i 


PN {a;^,q;zq’' ’') = 0 . 


It is now convenient to introduce the identities 










= (-ir i+E 




‘dS (A)' 


where 


Y[i^ + Pk x) = l + ^[bk (P) x'^] 

k^l k^l 

SO that the quantities bk (p) are defined as follows, 

s 

bi ii) = E , 




ki,k2 — l, ki^k2 
s 

h{P)= E {13k, Pk. f^kj , 

ki,k2:k3 — l, k\^k2, k2^k3, k^^ki, 


and so on, up to 


bs {§)= n ; 




(45) 


(46) 

(47a) 

(47b) 

(47c) 

(47d) 

(47e) 


12 












as well as the analogous identities 


r r r 



n (^“i)= n - 1 ) =n 

i=l k=l j=l 

( r '1 


where 

= (-1)’’ i-sy] I, 

r r 

(48) 


(1 -b Oj a;) = 1 + ^ [aj (a) x^] , 
i=i 1=1 

(49a) 

so that the quantities aj (a) 

are defined as follows. 



r 

ai (a) = (“a) > 

r 

(49b) 


(ii) 'y ^ (®ii 11 ^ 2 )) 

ii.i2=i. h/l2 

(49c) 


r 

03 {(£) = y ' (nji ’ 

ilj2j3 = l. j2^j3, 

(49d) 

and so on, up to 

T 



iy) = n (“ 1 ) • 

i=i 

(49e) 


Hence the ^-difference equation (HSl) satisfied by the generalized basic hypergeometric polynomial pn (z) (see 
(l3a|) and ([T])) reads as follows (see (I44bl) h 


^1+E 

k=l L 


(-f)' 


PN (z) 


-z {1-q ^S) (-!)’■ ® («) } PNizq^^ ^)= 0, 


(50a) 


hence 


{l - ^ + E [(-?)■" h {§) (<5" - I PN {z) 


- (-!)’■ ® z ll-q ^6+ (“ly («) - 9 > PAT (z "■) = 0 , (50b) 


hence 


Pn{z)-Pn{z q) + ^{-q) ^bk{^ [pN {z q'") - Pn {z q^^^)] 


fc=l 


- (-1)’'-^ z [pn {z - q-^PN {z 9*-’’+')] 

(-1)’'“'* z I ^ [(-1)^ Oj (a) PN [z 9*“’’+^) - q~^ ^ [(-1)^ Oj (a)PAT (z q^-'^+3+^) > = 0 . (50c) 




j=i 


Next, let us look at this formula—which is a polynomial equation in z of degree N + 1—at z = where 
is one of the N zeros of the polynomial pn (z), see (Hal) . It is then plain that there hold the N algebraic 
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equations 


-PW (Cn <?) + X! [PN {Cn q'") - PN {Cn 

k^l 

-{-ir‘ C„ [pA-ftn?—')-g-''p»(C„<i-'+‘)l 


f ^ 


(-!)'■■' c„ E 

(-1)^ Qj (a) PN (Cn 9" 

[j=l, jAr-s 



i=i, i/i'-s-i 

n=l,2,..,N. (51) 

And it is then immediately seen that the insertion in these formulas of the expression (|4bp of the monic 
polynomial pN {z) in terms of its zeros yields ([71). Proposition 2.1 is thereby proven. 
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